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Correspondence between frequency and energy spectra and biorthogonality conditions for the 
excitations of Bose-Einstein condensates described by Gross-Pitaevskii model have been derived 

■ selfconsistently revealing several novel aspects originating in nonselfadjoitness of the Bogoliubov 
' operator. It has been demonstrated that frequency resonances of the excitations with positive and 

negative energies can lead to their mutual annihilation and appearance of the collective modes with 
complex frequencies and zero energies. Conditions for the avoided crossing of energy levels have 
also been discussed. General theory has been verified both numerically and analytically in the weak 

■ interaction limit considering an example of vortices in a binary mixture of condensates. Growth of 
qq ' excitations with complex frequencies leads to spiraling of the unit and double vortices out of the 

condensate center to its periphery and to splitting of the double and higher order vortices to the 
unit ones. 
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I. INTRODUCTION 



o 

Recent observations of the quantized vortices in a dilute ultracold gas of 87 Kb atoms [|l]J^] are spectacular evidences 
of the superfluid properties of atomic gases below transition temperature for formation of Bose-Einstein condensate 
(BEC). First experimental results were obtained, using method suggested by Williams and Holland ||, in a binary 
mixture of the hyperfine states of 87 i?6 Q and more recent experiments Q have demonstrated vortices and vortex 
£j ■ lattices in an optically stirred single component condensate. 

Gross-Pitaevskii (GP) equation [Q is a widely used approximation to describe dynamics of the quantized vortices 
in the superfluid component of Bose gases at temperatures below critical. The GP equation falls into a general 
class of Hamiltonian nonlinear systems. The theory of stability of stationary solutions (equilibria) of such systems 
is well established at the moment, see e.g. ||. However, a certain gap between formal mathematical knowledge and 
' applications to physical examples still exists and BEC field is not an exception. One of the controversial examples 
in the BEC context is complex or imaginary eigenfrequencies in the spectrum of elementary excitations. Existence 
of these frequencies was first pointed out by Bogoliubov in his original work jfj|, for more recent references see, e.g. 
[fr| — - Unjustified negligence by the nonselfadjoitness of the Bogoliubov equations often lead to the association 
between frequencies and energies in a way standard for quantum mechanics based on selfadjoint operators, thus 
admitting possibility of complex energies in a conservative system. 

However, rigorously proved theorem allowing clear physical interpretation of complex frequencies is known from the 

S general theory of Hamiltonian systems j^] . It states that excitations with complex frequencies can appear only as a 
result of the resonance between two excitations with positive and negative energies (l3) . The roots of this theorem go 
back to 19th century |L4| and it has been used for interpretation of instabilities with complex eigenvalues in plasma 
physics p5| and fluid dynamics |16|j . Motivated by recent experiments on observation of vortices in binary mixture 
of condensates |l| we will demonstrate that these vortices can have complex frequencies in their spectrum, thereby 
giving good practical ground for selfconsistent theoretical interpretation of this phenomenon in the BEC context. 

Properties of vortices in the single component magnetically trapped ultracold gases have been subject to the 
intensive theoretical investigations, see e.g. [9,17-32], which significantly extended classical wor ks |JH dealing with 
spatially unbounded case. The properties of unit vortices in the two component condensates have also been studied, 
in parallel and independently from this work, by Garcia- Ripoll and Perez-Garcia |l(]Jll| ■ Richness of the dynamics of 
the two-condensate system and different approaches to the problem have led only to few overlaps which are outlined 
where appropriate. 

In the next section we introduce coupled GP equations and briefly describe their general properties. Then, in section 
III, we derive Bogoliubov equations for excitations, clearly specifying differences between frequency and energy spectra 
of the excitations. We also explain scenario of appearance of the excitations with complex frequencies and show that 
they have zero energies. In section IV we verify validity of general results presented in Sec. Ill using perturbation 



*URL: www://http. cnqo.strath.ac.uk/~dmitry 



1 



theory in the limit of weak interaction and direct numerical study of Bogoliubov equations. In Sections IV, V we 
describe how long term dynamics of unit and higher order vortices can be interpreted using linear Bogoliubov theory. 



II. GROSS-PITAEVSKII EQUATIONS 

Studies of superfluid mixtures using coupled GP equations have long history and attracted significant recent activ- 
ities, see e.g. ju[r|MH,0,Bi] [58| and references therein. Following these works we assume that wave functions V'1,2 of 
two-species condensate inside an axial harmonic trap obey equations 

h 2 1 
ihdt^i = -— V 2 V;i + -mfl 2 (r 2 + erVUi 
2m 2 

+ (u U |Vl| 2 + Ul 2 |V2| 2 )V'l, (1) 

h 2 - i 

iHdtih = -^V 2 V 2 + -mfl 2 {r 2 + a 2 z 2 )ip 2 
2m 2 



+ (m22|V j 2| 2 + M2i|^i| 2 )V': 



2- 



where for simplicity we have neglected by possible differences of atomic masses m 1,2 = m and trap frequencies, 
^1,2 = fJj ci,2 — f, V = i x d x + iydy + i z d Zl r 2 — x 2 + y 2 . Coefficients Uij — Airh 2 ay / m characterise intra- and 
inter-species interaction with corresponding two-body scattering lengths an ^ a 22 and 012 = & 2 i- 

At this point we introduce dimensionless time and space variables t = £lt and [x, y, z) = (x,y,z)/a,ho and normal- 
isation for the wave functions ^12 = a^ 2 ?/;i, 2 /v^Vi, where a^ = \Jtij (mil) is the harmonic oscillator strength and 
iVi^a = J dV\ipi.2\ 2 are the numbers of particles. We will consider the quasi-2D model to simplify our numerical study. 
This approximation was previously used in several works, see e.g. p|-pj|, and it is applicable not only for pancake 
traps, (7» 1, but also captures main qualitative features of spherical traps. To make further reduction of Eqs. (|l|) 
we redefine the wave functions once more: 
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1/4 *i,2(i,y,i)e- CTi;2 / 4 e- Vl - 2 *- l,T * /2 , (2) 



where /Lti 2 are the chemical potentials. Dropping tilde we find that equations for $ 12 and 2 can be put into 
Hamiltonian form 

m 

id4 + fj— = 0, (3) 
Sip* 
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H= J dxdy(\Wi\ 

+ (1>-Mi)l*i| 2 - 
9, 



|W 2 | 2 



(V- 1*2)1*2] 



(4) 



-|[/9n|*i| 4 + /3 22 |*2| 4 + 2/3 12 |^ 1 | 2 |* 2 | 2 ]). 

where H is the energy functional (or Hamiltonian), V = i x d x + i y d y , V = r 2 /4 is the harmonic potential, g is the 
interaction parameter g = E^/naNiau/ahoj P12 = 012/011, and /3 22 = a 22 /an. (3\\ = 1 and it is left in the equations 
for the sake of the symmetry. 

Invariancies of H with respect to the infinitesimal rotations and two parameter phase transformation 



(*i,* 2 )->(tfie* 1 ,*2e < * , ) 1 

result in conservation of the total angular momentum and of the total number of particles in each component. 
Radially symmetric stationary states of the condensate (equilibria) can be presented in the form 



(5) 



* J = A j (r)e iL ' e ) j = 1,2 



(6) 
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where 6 is the polar angle and Aj are real functions. Using method suggested in only states with vortex in one 
of the components can be created and therefore we will consider below only cases with L 2 > and L± = 0. Functions 
Aj{r) were found numerically using Newton method. Chemical potentials /ii.2 were found from the normalization 
conditions 

2tt J rdrA\ = 1, 2tt J rdrA\ = ^ = N. (7) 



III. FREQUENCY AND ENERGY SPECTRA OF COLLECTIVE EXCITATIONS 



A. Bogoliubov equations and frequency spectrum 

To study spectrum of BEC at equilibrium we linearise Eqs. (0) using substitutions 



(A,- (r) + &(»•, 0,t))e 



iL j 



(8) 



where fj are small and complex. Assuming that excitation are periodic in 9 with period 2n we expand fj into Fourier 



scries: 



no \ 



(9) 



where I = 0,±1,±2, . . . Then 



idtWi+fjHiW^O 



is the set of linear partial differential equations resulting from the substitution of 

Wi = {U ll ,V a ,U l 2,V l2 ) T , 



into Eq. 



(10) 
Here 



Hi = 



A,i gPnAf 9P12A1A2 gf3i 2 AiA 2 
gP 1 2A 1 A 2 gPi2A!A 2 L1.2 9P21AI 



gPi 2 AiA 2 gPi2A x A 2 g$<iiA\ 



1:2 



is a self-adjoint operator and 



£1,3 



1 d d 



r dr dr 

\2 , „a 42 



+ ^(Lj+l) 2 + V - N 



+2gP jj Aj + gpuAjr, 3 = 1, 2; f - 2, 1. 

Frequency spectra of rjHi are discrete, providing V =/= 0, therefore phonons, strictly speaking, are absent in the trap 
geometry. In accord with standard terminology [^j39|, all spatially bounded elementary excitations can be called 
collective excitations (or collective modes) of an equilibrium under consideration. Linearised equations for excitations 
in a Bose gas, similar to Eq. (|l(]), were first derived by Bogoliubov and expansion (^) was first applied in the context 
of the vortex excitations by Pitaevskii To find frequency spectrum and collective modes we need to solve set of 
the eigenvalue problems 



(11) 



ffHi are non-self-adjoint operators and therefore complex frequencies are not forbidden. If wi n is an eigenvector of 
ffHi with eigenvalue loi„ it is selfevident that wf n is also an eigenvector with eigenvalue ixi\ n and it can be shown that 



ffH-i has eigenvectors W-i n = fwi n and w*_ ln with eigenvalues —uoi n and 
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respectively. Here 
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Thus spectrum of ffrt_i can be obtained by reflection of the spectrum of ffhti with respect to the line Reu> = in 
the plane (Reu>,ImLj). In other words, it means that purely real or purely imaginary frequencies of the elementary 
excitation exist in pairs and complex frequencies exist in quartets and that 

Tr{f){Hi + H-i)} = ^{Ujln + LJ-ln) = 0. (12) 
n 

Any equilibrium state of the condensate is spectrally stable || if its spectrum is real. If there is at least one frequency 
with negative imaginary part then corresponding collective mode will grow in time destabilising the equilibrium, which 
is called spectrally unstable. 



B. Biorthogonality 

Eigenmodes of fffii are biorthogonal to the modes of the adjoint eigenvalue problem p0| , i.e. 

(win,a ln ,) = 0, (13) 

where n ^ n' and <z;„ obey 

{fjHlfdin = Hlfjain = UJi n ai n , (14) 

and (6, c) = 2tt ^2 k J °° rdrb* k Ck for any b, c. Factor 2tt is introduced to mimic integration over 9. The key feature 
of our model, originating in its Hamiltonian structure, is that transformation linking wi n and its adjoint a;„ can be 
found in explicit and simple form. If wi n is an eigenmode of ffHi with frequency loi„ then it can be checked that 
f)w* n and fjwin are eigenmodes of Jiif] with eigenvalues, respectively, w ; * n and oji n . The mode adjoint to wi n is fjwf n , 
therefore if Imuji n ^ then biorthogonality condition ( |l3| ) implies 

(wi n> fr&ln) = 0- (15) 

If wi n is real, then wi n is also real and {win,f]win) ^ 0. Normalisation constant can always be chosen in such a way 
that 

\{wi n ,rfuJi n )\ = 2. (16) 



Convenience of making the left hand side of Eq. (|6|) equal to 2 will become clear below, when Eqs. (19), (20) 
for energies of elementary excitations are derived. Eq. (1^) makes it explicit that inner product (wi n , fjwi n ) can be 
either positive or negative. This point often remains silent if one derives conditions similar to (nq) as part of the 
diagonalisation procedure of the second-quantised Hamiltonian disregarding eigenmodes with negative and zero values 
of (winjfpjjin) |||. The fact that inner product of a real eigenmode with its adjoint can be negative, is different from 
standard quantum mechanics based on the selfadjoint operators. These difference can have series of consequences and 
one of them is that energy levels are not necessarily linked to the eigenfrequencies according to the standard rule hw, 
see below. Note, that origin of the nonsclfadjoitness in our case is the nonlinearity of GP equations. If particles in 
the condensate do not interact, g — 0, then fjTLi become diagonal and self-adjoint. 



C. Energy spectrum 

As a prelude to calculation of energies of the elementary excitations it is instructive to introduce notion of the 
nonlinear stability, i.e. stability under the full nonlinear dynamics. According to Dirichlet's theorem nonlinear 
stability is ensured if an equilibrium state under consideration is either minimum or maximum of the functional 
H. Note here, that excitations change number of particles in the equilibrium state, therefore it was convenient to 
introduce energy functional H, which is actually the so-called modified energy |33| , i.e. it is the energy functional for 
Eq. (1) modified by addition of the number of particles integrals, see terms proportional to fij in Eq. (§). Let us 
assume that wi{r)e lW is a small initial perturbation of an equilibrium state of the condensate, then 

H = H + ^{w l ,n l w l ) + ..., (17) 

where H is the energy of the perturbed equilibrium and Hq is calculated at the exact equilibrium. The equilibrium 
is nonlinearly stable if eigenvalue problems 
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(18) 



have all their eigenvalues cither negative or positive, except zero eigenvalues generated by continuous symmetries. 
Spectral instability implies nonlinear one, and nonlinear stability implies spectral one, but not vise versa ]^|. In the 
nonrotating traps all the higher order states of the condensate including vortices are not the local extrema of the 
energy p9| , p3p5|j27[ | . Let us stress, however, that their nonlinear instability can not be guaranteed by this fact alone 
and requires separate consideration. 
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FIG. 1. Numerically calculated frequency (a) and energy (b) spectra corresponding to the collective modes with / = ±1 of 
the spectrally stable (left panel) and spectraly unstable (right panel) unit vortices: Lj = 0, L% = 1, g = 250, N = 1, other 
parameters for the left/right panel correspond to the vortex in the state 1/state 2, see text after Eq. (32). The correspondence 
between frequency and energy spectra is obvious, see Eqs. (^,^l|). (c) Spectrum of the eigenvalue problems Ti±i/3„ = a n 0n', 
parameters as for (a),(b). 



If wi in Eq. ( |l7| ) is an eigenmode u>/„ of fffii then 

1 



ejr, 



(wi n ,HiW ln ) = -Ul n {wi n ,flWi n ) 



2 „w 2 

measures energy of this mode. Assuming that cj;„ is real and using biorthogonality conditions (|l^) one gets 

(In = LUl n sign(wi n , fjWln). 



If Imujin ^ then Eq. <ha) implies 



Energy in physical units is given by ei n h£l. 
It is readily demonstrated that 



{Wln,T]Wln) = ~ {TW ln , T]TW[ n ) = -{W-l n , T]W-l n ) 



(19) 

(20) 
(21) 

(22) 



therefore modes with frequencies u>i n and u)-i n = —u>i n have their energies equal in sign and value. 

Eigenfunctions /3 m ; are orthogonal (AmjAm') = <$mm' and form complete basis. Therefore wi n in Eq. (|l^ ) can 
be expanded in terms of /3; TO , which gives e;„ = J2 m a im(wi n , Pi m } 2 ■ Thus an equilibrium can have collective modes 
carrying energy with opposite signs only if it is not a local extremum of the energy functional. 
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D. Resonances and instabilities, crossings and avoided crossings 



Instabilities with imaginary eigenvalues have been found to dominate dynamics of homogeneous BEC with attractive 
interaction |(| and mixture of two condensates with repulsive interaction • Higher order vortices and dark solitons 
in single component trapped condensates can be unstable with respect to modes with complex frequencies [^,^2[- In 
our problem instabilities with complex eigenvalues are the most important ones also. Therefore it is desirable to have 
criterion or theorem allowing interpretation and prediction of these instabilities. Such theorem is actually known from 
the general theory of the Hamiltonian dynamical systems [^) and, in our context, it can be reformulated as: (i) Sign 
of the energy of a collective excitation preserves as parameters vary as long as there is no frequency resonance with 
another excitation, (ii) Condensate at equilibrium can lose spectral stability as parameters vary only in two ways: 
either via frequency resonance of two elementary excitations with positive and negative energies or by resonance at 
zero frequency. Proof of these results || is based on the fact, that transition from spectral stability to instability can 
not violate energy conservation law: dfH = 0. 

From this theorem and preceding considerations one can conclude that complex eigenvalues in the spectrum of the 
vortices in Bose-Einstein condensate can appear only due to mutual annihilation of the collective excitations with 
positive and negative energies. The latter can coexist only for nonlincarly unstable equilibria. The theorem does 
not forbid for frequencies of two excitations with either the same or opposite signs of the energies simply cross each 
other without change of the spectral stability. In our example it typically happens when corresponding eigenmodes 
remain orthogonal at the exact resonance. If, however, eigenmodes of the colliding excitations start to compete for the 
same direction in the functional phase space and become degenerate at the resonance then crossing is not a generic 
scenario. If the signs of the energies of the excitations are opposite, then quartet of complex frequencies appears 
upon passing the resonance. Alternatively, if the signs are the same, the exact resonance can not be achieved and it 
becomes replaced by the so-called avoided crossing of the energy levels §. 

Figs. 1(a), (b) show numerically calculated frequency and energy spectra of the collective excitations with I = ±1 
for spectrally stable (left panel) and spectrally unstable (right panel) unit vortices (Li = 0, L\ — 1). The negative 
energy excitation is clearly seen in the spectrally stable situation. At the transition threshold to spectral instability 
this excitation and another one, having energy with the same absolute value but opposite sign, annihilate each other. 
This transition is accompanied by appearance of the zero energy excitations. Examples of frequency and energy 
evolution under the parameter variation resulting in instabilities, crossings and avoided crossing are given in the 
section IV. Figs. 1(c) show spectra of H±i illustrating that there are no qualitative changes in these spectra after 
appearance of complex frequencies. 

Possibility of observation of the collective modes with complex frequencies has caused some concerns and discussions 
[pL2| , p7|| . However, initial perturbations of the equilibrium state having nonzero projections on the adjoint mode fjwf n 
with Imuj* n > will lead to the ultimate growth of the corresponding mode wi n because it has Imcuin < 0. The 
consequences of this growth is in any way diminished by the fact that the energy of this mode is zero. Wealth of 
references and examples of instabilities with complex eigenvalues existing in other physical contexts can be found in 



E. Goldstone and dipole modes 

Infinitesimal variations of <j)j, see Eq. (||), generate two zero-energy eigenmodes (Goldstone modes) (Ai, — A\, 0, 0) T 
and (0, 0, A2, — A2) 7, belonging to the null-eigenspace of ffH.^. Harmonic trapping modifies spectrum in such a way 
that ffH.±\ acquire couple of parameter independent eigenvalues u> = ±1 with eigenfunctions 



w±u = 



dA 1 
dr 


T^iiAi) 


+ 2^1 


dAi 
dr 


±lL 1 A 1 ) 


-\rA x 


dM 
dr 


T l r L 2 A 2 ) 


+ \rA 2 


dA 2 
dr 


± \L 2 A 2 ) 


- ¥A 2 _ 



(23) 



where u> can take values ±1 for both eigenmodes. w±id are often called dipole modes and e±id = 1. Eqs. ( |23| ) 
generalise expressions previously derived for single-species condensates Existence of the dipole modes can also 
be associated with the Kohn's theorem flTdL[42|] . 
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IV. COLLECTIVE EXCITATIONS OF UNIT VORTICES: Li = 0, Li = 1 

We begin our analysis considering weakly interacting condensate: g -C 1, N ~ 1. In this limit potential energy due 
to harmonic trapping V strongly dominates over the interaction energy, which allows to make explicit calculations of 
the frequency and energy spectra of the elementary excitations. Calculations of the excitation spectra of the vortices 
and dark solitons in the singly component weakly interacting condensates have been recently done by several groups 
of authors [|l2|]25| l . However, these calculations lack analysis of the energy sign in a sense explained in the preceding 
section. 

We substitute asymptotic expansions /z^a = a4°2 + 9l*i\ + 0(g 2 ), ^1,2 = 4°2 + gA^\ + 0{g 2 ) into the stationary 
(dt = 0) version of Eqs. ^ and derive recurrent system of linear equations. In the zero approximation we have two 
uncoupled harmonic oscillator problems with eigenmodes 

4 0) = A (0) = M* e -^/*. (24) 

V27T V 7T 2 

Using solvability condition of the first order problem we find asymptotic expressions for the chemical potentials: 
Mi = 1 + 3(2/3ii + nA 2 )/(87r) + 0(g 2 ), M2 = 2 + g(n(3 22 + 12 )/{8n) + 0(g 2 ). 

Then we expand operators Hi, eigenmodes wi and frequencies u>i into the series Hi — Hi + gHi + 0(g 2 ), 

wi = w[ ^ + gw^ + 0(g 2 ), uji = lo^ + gco^ + 0(g 2 ). After substitution into Eq. ( |TT| ) in the first order we find 
standard equation 

(^P>- W p>)^ = (a,«-^)^. (25) 

ffftf 1 is diagonal and self-adjoint and all its eigenmodes and eigenvalues can be found explicitly. Then using solvability 
condition for Eq. ( p5[ ) one can find corrections for all frequencies and coefficients in the linear superposition of the zero 
approximation eigenmodes. Computer algebra makes technical realization of this plan by a straightforward exercise. 
We will present and analyse explicit analytical results only for the operator ffH-i considering vicinity of the spectral 
point u) = 1, because it contains information about origin of spectral instabilities of unit vortices. Equivalent analysis 
of r)Hi near uj — — 1 gives the same results. 

ffH^\ has three unit frequencies, = 1. Corresponding eigenmodes are 

6 1 = (0,l,0,0) T -^e-'' 2 /4 j 

6 2 = (0,0,l,0) T 4=e~ r2/4 , (26 ) 

v Zir 



&3 = (0,0,0,l) T -^e-'' 2 / 4 . 

Solvability conditions for Eq. ( ^5| ) lead to the characteristic determinant of the three by three matrix. We find that 
one of the three frequencies is associated with dipole mode w-id and that the other two are 

= 1 + 327 ( ~ 3/312 ~ NP22 ± ^)+ °^ 2 ), (27) 
R = (3/3i 2 + N0 22 ) 2 - 80 2 2 (N + 1). 

Corresponding unnormalised eigenmodes are 

wtt = V2N (n0 22 ~ fa ± VM) h 

+V2 {iN0 22 + [1 - AN]0 12 T v 7 !?) b 2 

+4A(/3i 2 - Pn) h + 0(g). (28) 
Pair of frequencies ( p7| ) becomes complex if R < 0, signaling spectral instability of the unit vortex. If R > 0, then 

(w ± 1 ,f t w ± 1 ) = (29) 
4 \R(N - 1) ± y/R{(l - 5N)0 12 + N(3 + N)0 22 } 
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and one can check that 



e± x = ±u% (30) 

i.e. modes with frequencies u>_\ and uZ\ have, respectively, positive and negative energies, see Fig. 2(c). One can 
also verify biorthogonality (ur^i, fjiffti) = = e i z l for R < 0, see Eq. (|T5|), and (w^ 1 , fjw^i) = for R > 0, see Eq. 
©• 

Assuming > and rewriting instability condition R < in the form 

§^<l(2v/2(iV+l)-3), (31) 

Pl2 N 

one can see that unit vortices are more stable if intraspecies interaction of atoms in the vortex containing part of the 
condensate is somewhat larger than interspecies interaction. The choices of scattering lengths corresponding to the 
experiment @ are: fa = 1.03 /0.97, fa = 1/0.97 (fa > fa) and fa = 0.97/1.03, fa = 1/1.03 (fa < fa). The 
former case corresponds to the vortex in the spin state {F = 1, m/ = —1} of 87 Rb and the latter to the vortex in the 
state {2, 2}. These two states will be called - state 1 and state 2. It is clear that for N = 1 Eq. ( pl| ) predicts instability 
for the vortex in the state 2 and stability for the vortex in the state 1, which supports results of the experimental 
observations [Q. 

Fig. 2 shows frequency resonance accompanied by the simultaneous mutual annihilation of excitations with positive 
and negative energies happening at some critical value of /?i 2 , see Eq. (|3l"|). In fact it models transition from the 
situation with vortex in the state 1 to the case with vortex in the state 2. Performing numerical studies for wide 
range of parameters, outside the validity region of analytical considerations, we have not been able to find regions of 
spectral instability of the unit vortex in state 1. Contrary, existence of instabilities of the vortex in state 2 due to 
exactly the same scenario, which is predicted in the weak interaction limit, can be readily demonstrated, see Fig. 3. 

It is important to stress, that if two condensates are decoupled, fiyi = 0, then the negative energy mode wZi belongs 
to the vortex containing component and the positive energy mode w_\ belongs to the vortex free component. Because 
of this separation instability is not possible for any values of g, which agrees with spectral stability of unit vortices 
in the singly component case reported in |J. Thus we can conclude that the vortex instability in our example has 
essentially two-component nature and its analog may also exist in the case when second component is a non-condensate 
one. 

Criterion similar to ([3l]), but without corresponding energy analysis, has also been independently obtained in p"l| ] 
using two mode approach, i.e. condensate wave functions \?i and ^ 2 have been presented as linear superposition of 
modes (Q) with time dependent coefficients and GP equations have been reduced to the set of ordinary differential 
equations for these coefficients jn|. However, this method fails taking into account an eigenmode proportional to r 2 , 
see 63 in (]2^), which makes an important contribution to the expressions for frequencies. Therefore Eqs. (|27j) and 
( |3l| ) are different from the corresponding results presented in JTl| ] . 

To illustrate crossings and avoided crossings in the spectrum of unit vortices we show in Fig. 4 positive part of the 
frequency spectra of f]TL±i and ffH±2- One can see numerous points, which on the first glance can be interpreted as 
crossings in the frequency spectrum. However, under the close investigation those of them which are marked by the 
open circles, turns to be the avoided crossing of the excitations with equal energy signs, see Fig. 5. 
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FIG. 2. (a),(b) Frequencies uCj vs /?i2: N = 2, g = 0.1, ,022 = 0.97/1.03. (c) Inner products {w_ 1 ,rjw_ 1 ) characterising 



energy signs of the corresponding excitations, sec Eq. (30). Doted lines are numerical solutions of fjlK-iWi 
lines correspond to Eqs. (28) and (30). 
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FIG. 3. Instability growth rate of the unit vortex in state 2 vs interaction parameter g. Full line - N = 1; Dashed line - 
N = 0.3; Dotted line - N = 9. 
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FIG. 4. Positive part of the frequency spectrum of the collective excitations with |/| = 1 (full lines) and \l\ = 2 (doted lines) 
of the unit vortex in state 1 vs interaction parameter g: N — 1. Open circles mark the avoided crossings (i.e. avoided frequency 
resonances) of the excitation with positive energies. 



FIG. 5. Dependencies of energies of two elementary excitations with I — 2 vs g showing details of the rightmost from the 
avoided crossings marked in Fig. 4. 



The unstable modes of the vortex in the state 2 are of the dipole type, i.e. |Z| = 1, therefore their growth leads to the 
initial displacement of the vortex from the trap center. Displaced vortex carries less of the total angular momentum, 
compare to the momentum of the vortex positioned at the trap center. Lack of the angular momentum is compensated 
by two factors. First, vortex acquires a nonzero tangential velocity and therefore its trajectory is actually a spiral, 
which is similar to the dynamics of an optical vortex displaced from the center of a gaussian beam fj^j] . Second, 
angular momentum, and vortex itself, become gradually transferred into the second condensate component, which 
was first demonstrated in |^0| , p| . 

If drift instability is absent, then dissipative effects still can result in the vortex drift, which was predicted for the 
single component condensates by several authors [ ^9|j32| . Presence of the second condensate opens a channel for the 
energy and angular momentum transfer from the vortex into the vortex free component. It can be seen that vector 
&i, being excited by the instability, provides a channel for this transfer. Thus drift instability can also be interpreted 
as due to dissipation of the energy and momentum by the vortex free condensate component. 

In the limits N 3> 1 and N <C 1 our model can be approximately considered as a single component condensate with 
(N ^> 1) or without (N <C 1) vortex. Unit vortex and ground state of the single component condensate are known 
to be spectrally stable. Therefore drift instability disappears in both limits, see Fig. 3. Increase of g for fixed N also 
results in the suppression of the instability, see Fig. 3, which means that not only relative, but also absolute increase 
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V. DRIFT OF UNIT VORTICES 
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of the number of atoms in the vortex free component stabilizes the condensate. 

VI. DRIFT AND SPLITTING OF HIGHER ORDER VORTICES 

Considering higher order vortices one can expect to find instability scenario resulting in their splitting into unit 
vortices. This scenarion is expected to be due to growth of the collective modes with \l\ > 1. However, as we will see 
below the drift instability linked to the dipole like modes, \l\ = 1, also can be presented. It leads to displacement of 
the whole vortex from the trap center without splitting, at least at the onset of the instability. 

Both drift and splitting instabilities of the higher-order vortices appear as a result of the frequency resonances of 
the elementary excitations with negative and positive energies, similar to the case of unit vortices. N <C 1 corresponds 
to the vortex free condensate and therefore both instabilities disappear in this limit. In the limit N ^> 1 only drift 
instability is suppressed and one can recover periodic in N bands of the instabilities with complex frequencies and 
\l\ > 1 similar to the results reported for higher order vortices in a singly component condensate Note, that 
splitting itself was not explicitly demonstrated in ||. It is also interesting to note that higher-order vortices in the 
free, V — 0, single component condensate are spectrally stable [j45[ . Thus splitting can be considered as induced by 
the trapping. The vortex free condensate component plays crucial role in the drift instability, but will the latter one 
be presented without trapping or not remains an open problem. 
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FIG. 6. Growth rate of the drift |/| = 1 (full lines) and splitting \l\ — 2 (dotted lines) instabilities of the double vortex in the 
state 2 vs N for g = 900. 

Double vortices: L\ = 0, L2 = 2 

Considering double vortex we have found that it can be unstable with respect to the |Z| = 1,2 excitations. The 
vortex in the state 1 has been found surprisingly stable. One has to take relatively small values of g and large N to find 
splitting instability. Vortex in the state 2 is more unstable in a sense that splitting exists already for N ~ 1, see Fig. 

6. As it is evident from Fig. 6 either drift or splitting instability can dominate vortex dynamics. If drift instability 
is dominant, then vortex first gets displaced from the trap center and only then splits into the unit ones, see Fig. 

7. The latter happens due to the curved background which breaks cylindrical symmetry with respect to the vortex 
axis. After the splitting vortices remain close to each other and move towards the condensate periphery. Results of 
the numerical simulation of CP equations (1) presented in this Section were obtained starting from equilibrium states 
perturbed by random noise with amplitude ~ O.ObAi^- 

The dynamics is quite different when splitting instability is dominant, see Fig. 8. In this case unit vortices appear 
straight at the onset of the instability development and spiral out of the condensate center being always positioned 
symmetrically with respect to it. After a certain period of time vortices move back to the trap center and condensate 
state close to the initial one is restored, see Fig. 8, then the cycle is repeated with gradually worsening degree of 
periodicity. 

During development of the instability angular momentum and vorticity become partially transferred into the second 
condensate. Analysis of the transverse profiles of the phases corresponding to the density profiles shown in Fig. 8, has 
revealed that black spots appearing in the second condensate are indeed unit vortices, not the density holes without 
topological structure. 
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FIG. 7. Development of the drift instability of the double vortex in the state 2: g = 900, N = 2.5. Top - \^i(x,y)\ 2 , bottom 
- 1*2(2;, y)\ 2 ■ Time interval between snapshots is 20. First snapshot corresponds to t = 60. 
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FIG. 8. Development of splitting instability of the double vortex in the state 2: g = 900, N = 0.7. Top - |*i(a;, y)\ 2 , bottom 
- 1*2(2;, y)\ 2 ■ Time interval between snapshots is 20. First snapshot corresponds to t = 30. 

Triple vortices: L\ = 0, 1/2=3 

Reach variety of beautiful vortex lattices can be found considering instabilities of vortices of the order 3 and higher. 
This reachness can be understood in terms of spatial profiles of the unstable collective modes. E.g. triple vortex has 
been found unstable with respect to the perturbations with \l\ = 1,2,3. All components of the \l\ = 2 excitations 
are equal to zero at the trap center. Therefore one can expect that growth of this mode will develop into a spatial 
structure preserving vortex at the trap center. Contrary some of the components of the \l\ = 3 modes have humps at 
the center. Therefore their growth should repel all unit vortices out of the center, which leads to the breaking of the 
triple vortex into a triangular structure of the unit vortices moving away from the trap center. Both instabilities can 
be found for the same parameters values and can have close growth rates. Therefore the winning mode is selected 
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through the process of complex competition, see Figs. 9,10. 
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FIG. 9. |Zj = 2 instability of the triple vortex in the state 2: g = 900, N = 0.9, t = 70. Left - |*i (x, y) \ 2 , right - \^ 2 (x,y) 
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FIG. 10. \l\ = 3 instability of the triple vortex in the state 2: g = 900, N = 0.9, t = 85. Left - |*i(:r,y)| 2 , right - \^ 2 (x,y)\ 2 . 

VII. SUMMARY 

We have described general ap proach to stability of equilibrium in BEC using Bogoliubov theory and GP equation. 
Biorthogonality conditions (|lq), ( |l6|) and correspondence between frequency and energy spectra of the elementary 
excitations (pfl|), ( pi] ) have been derived selfconsistently from first principals revealing several novel and conceptually 
important aspects originating in nonseldjoitness of the Bogoliubov operator. 

It has been demonstrated that frequency resonances of the excitations with positive and negative energies can 
lead to their mutual annihilation and appearance of the collective modes with complex frequencies and zero energies. 
Conditions for the avoided crossing of energy levels have also been discussed. General theory has been verified both 
numerically and analytically in the weak interaction limit considering an example of vortices in a binary mixture of 
condensates. 

Growth of excitations with complex frequencies leads to the two main scenarios of the instability development. 
First one is the spiraling of unit and double vortices out of the condensate center to its periphery. Second scenario is 
the splitting of the double and higher order vortices into unit ones. Absolute and/or relative increase of the number 
of particles in the vortex free condensate component have been found to have stabilizing effect. 

ACKNOWLEDGMENTS 

Author is particularly grateful to referees for their critical and very helpful comments. He also acknowledges 
discussions with S.M. Barnett and W.J. Firth. Numerical part of the work was significantly speeded up due to 
access to the computer equipment obtained via U.K. EPSRC grant GR/M31880 and assistance of G. Harkness and 
R. Martin. 



[1] M.R. Matthews, B.P. Anderson, P.C. Haljan, D.S. Hall, C.E. Wieman, and E.A. Cornell, Phys. Rev. Lett. 83, 2498 (1999). 



13 



K.W. Madison, F. Chevy, W. Wohlleben, and J. Dalibard, Phys. Rev. Lett. 84, 806 (2000). 
J.E. Williams and M.J. Holland, Nature 401, 568 (1999). 

L.P. Pitaevskii, Zh. Exp. Teor. Fiz. 40, 646 (1961) [Sov. Phys. JETP 13, 451 (1961)]. 

R.S. MacKay, Stability of equilibria of Hamiltonian systems, in 'Hamiltonian Dynamical Systems', eds. R.S. MacKay and 
J.D. Meiss (Hilger,Bristol,1987), page 137. 
N. Bogolubov, J. Phys. 11, 23 (1947). 

E.A. Goldstein and P. Meystre, Phys. Rev. A 55, 2935 (1997). 
E. Timmermans, Phys. Rev. Lett. 81, 5718 (1998). 

H. Pu, C.K. Law, J.H. Eberly, and N.P. Bigelow, Phys. Rev. A 59, 1533 (1999). 
J.J. Garcia- Ripoll and V.M. Perez-Garsia, Phys. Rev. Lett. 84, 4264 (2000). 
J.J. Garcia-Ripoll and V.M. Perez-Garsia, Phys. Rev. A 62, 033601 (2000). 



D.L. Feder, M.S. Pindoza, L.A. Collins, B.I. Schneider, and C.W. Clark, |sond-mat/0004449 . 
In the mathematical literature the energy sign is called signature of the eigenvalue or Krein sign. 

K. Weierstrass, in Weierstrass: Mathematishe Werke, vol. 1 (Berlin, 1894); reprinted (Hildesheim and New York: Olms 
and Johnson, 1967) page 233. 

A. Hasegawa, Plasma Instabilities and Nonlinear Effects (Springer, Berlin, 1975), page 18. 
R.S MacKay and P.G. Saffman, Proc. R. Soc. A 406, 115 (1986). 

M. Edwards, R.J. Dodd, C.W. Clark, P.A. Ruprecht, and K. Burnett, Phys. Rev. A 53, R1950 (1996). 
F. Dalfovo and S. Stringari, Phys. Rev. A 53, 2477 (1996). 

D. S. Rokhsar, Phys. Rev. Lett. 79, 2164 (1997). 
S. Sinha, Phys. Rev. A 55, 4325 (1997). 

R.J. Dodd, K. Burnett, M. Edwards, and C.W. Clark, Phys. Rev. A 56, 587 (1997). 
F. Zambelli and S. Stringari, Phys. Rev. Lett. 81, 1754 (1998). 

A. A. Svidzinsky and A.L. Fetter, Phys. Rev. A 58, 3168 (1998); A. A. Svidzinsky and A.L. Fetter, Phys. Rev. Lett. 84, 
5919 (2000). 

E. Lundh, C.J. Pethick, and H. Smith, Phys. Rev. A 58, 4816 (1998). 
M. Linn and A.L. Fetter, Phys. Rev. A 60, 4910 (1999). 

D.L. Feder, C.W. Clark, and B. Schneider, Phys. Rev. Lett. 82, 4956 (1999). 
J.J. Garcia-Ripoll and V.M. Perez-Garsia, Phys. Rev. A 60, 4864 (1999). 
D.L. Feder, C.W. Clark, and B. Schneider, Phys. Rev. A 61, 011601 (2000). 

B. Jackson, J.F. McCann, and C.S. Adams, Phys. Rev. A 61,013604 (2000). 
T. Isoshima and K. Machida, Phys. Rev. A 59, 2203 (1999). 

S. Stringari, Phys. Rev. Lett. 82, 4371 (1999). 

P.O. Fedichev and G.V. Shlyapnikov, Phys. Rev. A 60, R1779 (1999). 

A. L. Fetter, Ann. Phys. 70, 67 (1972). 

W.B. Colson, A.L. Fetter, J. Low Temp. Phys. 33, 231 (1978). 
T.-L. Ho and V.B. Shenoy, Phys. Rev. Lett. 77, 3276 (1996). 

C. K. Law, H. Pu, N.P. Bigelow, and J.H. Eberly, Phys. Rev. Lett. 79, 3105 (1997). 

B. D. Esry and C.H. Greene, Phys. Rev. A 59, 1457 (1999). 
P. Ohberg, Phys. Rev. A 59, 634 (1999). 

P. Nozieres and D. Pines, The Theory of Quantum Liquids. Vol. II. Superfiuid Bose Liquid (Perseus Books, Cambridge, 
1999), page 169. 

P.M. Morse and H. Feshbach, Methods of Theoretical Physics, Part 1 (McGraw-Hill, NY, 1953), page 884. 

W. Kohn, Phys. Rev. 123, 1242 (1961). 

J.F. Dobson, Phys. Rev. Lett. 73, 2244 (1994). 

D. V. Skryabin, Physica D 139, 186 (2000); D.V. Skryabin, Phys. Rev. E 60, 7511 (1999); D.V. Skryabin and W.J. Firth, 
Phys. Rev. E 58, R1252 (1998); ibid., 3916 (1998). 

Y.S. Kivshar, J. Christou, V. Tikhonenko, B. Luther-Davies, and L.M. Pismen, Opt. Commun. 152, 198 (1998). 
I. Aranson and V. Steinberg, Phys. Rev. B 53, 75 (1996). 



14 



